d + zd'-wave Superconducting States in Graphene 
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We show that effective superconducting orders generally emerge at low energy in the supercon- 
ducting state of graphene with conventionally defined pairing symmetry . We study such a particular 
interesting example, the d x 2_ y 2 +id' xy spin singlet pairing superconducting state in graphene, which 
can be generated by electronic correlation as well as induced through a proximity effect with a 
d-wave superconductor. We find that effectively the d-wave state is a state with mixed s-wave 
and exotic p + ip-w&ve pairing orders at low energy. This remarkable property leads to distinctive 
superconducting gap functions and novel behavior of the Andreev conductance spectra. 

PACS numbers: 74.45.+c,74.78.Na 



Graphene is a single layer of hexagonally coordinated 
carbon atoms which has recently been isolated Due to 
its special lattice structure, the low energy part of its en- 
ergy spectrum is characterized by particle-hole symmet- 
ric linear dispersions around the corners of the hexag- 
onal Brilliouin zone (BZ). This band structure is re- 
sponsible for many new properties of this 'relativistic' 
condensed matter system, such as an abnormal quan- 
tum Hall effect 0, HQ, minimum conductance 4, 51, and 
possibly even an experimental realization of the Klein 
paradox 

Recently, a novel concept called specular 
Andreev reflection was proposed for a nor- 
mal/superconducting(N/S) graphene interface in 
the context of a conventional s-wave pairing super- 
conducting stateQ- Later, an unusual oscillation of 
the quantum conductance through an N/I/S junction 
was predicted [1, Q. The possible superconducting 
pairing orders have also been studied. In Ref. 1 10 . 
by including strong electronic correlations, the mean 
field search shows that d x i_ y i + id' -wave pairing 
symmetry is favored, similar to the superconducting 
state in the triangular lattice which is believed to be 
of d x 2_ y 2 + id' symmetry [Til In Ref. [H, an 

exotic p + ip-w&ve superconductor with spin singlet 
bond pairing was suggested at the mean field level 
and possible phonon- or plasma-mediated mechanisms 
were discussed. On the other hand, experimentally, 
superconducting states i n g raphene have been realized 
by proximity effect [3, Il5l [la ] through contact with 
superconducting electrodes. 

The peculiar physics in graphene is the unusual lin- 
ear and isotropic dispersion of the low energy excitations 
around the Dirac points. In this Letter, we show that be- 
cause of the existence of the Dirac points, conventional 
pairing order parameters can lead to the emergence of ex- 
otic pairing states in the low energy effective description. 
The p+ip superconducting order of Ref. j_13| is precisely 
such an example of an effective, low energy supercon- 
ducting order, arising in that case from a more conven- 
tional extended s-wave pairing. Here, we study a partic- 
ularly interesting superconducting state in graphene, the 



d x 2_ y 2 +id xy spin singlet pairing superconducting state, 
which can be generated by electronic correlation[l(| as 
well as induced through a proximity effect with a d-wave 
superconductor on top of or underneath the graphene 
layer. We find that the d-wave state is effectively a mixed 
s-wave and exotic p + ip-wave pairing state at low en- 
ergy. The mixture of both s-wave and p + ip wave leads 
to unique properties of the excitation spectrum and An- 
dreev conductance spectra. The excitation spectrum is 
gapless at half-filling and is gapped away from half filling. 
The gap is equal to the chemical potential near half fill- 
ing, and it saturates as the chemical potiential is moved 
above the energy scale set by pairing strength. The nor- 
malized Andreev conductance in the limit of zero bias 
voltage is a smooth function of the chemical potential, 
which starts from 2 at half filling and drops smoothly 
to 4/3 at large doping, unlike that in the s-wave pair- 
ing states where it almost remains at a constant value, 
4/3 (see Fig. E^c)). This is a signature of d x 2_ y 2 + id xy 
pairing in graphene. 
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FIG. 1: (Color online.) Phase (blue number) of singlet bond 
pairing function on the graphene lattice which preserves the 
translational and rotational symmetry of the honeycomb lat- 
tice and is d + id' type under point group D%. The red vectors 
<5 a (a=l,2,3) denote nearest-neighbor inter-sublattice connec- 
tions. 



General pairing symmetry in graphene: Although the 
crystal point group of graphene is Dqh , the pairing sym- 
metry of the superconducting orders in a two dimen- 
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sional graphene sheet is governed by Dq, which includes 
four onc-dimcnsional irreducible representations, A\ 2 
and Bi^, and two two-dimensional irreducible represen- 
tations, E\_2- Among these representations, the A\,E\ 
and E2 representations describe s-wave, p-wave and d- 
wave pairing symmetries, respectively. Therefore, the 
spin singlet s-wave and d-wave pairing are described by 
the Ai and E2 irreducible representations. We can un- 
derstand the pairing symmetry further by considering the 
exchange symmetry between the A and B sublattices. 
The D 6 group is a direct product of its two subgroups 
Cs v and Z 2 , i.e. D 6 = C% v &)Z 2 , where Z2 describes the 
exchange operations between the A and B sublattices. 
The A\ and E2 representations of Dq are symmetric un- 
der exchange of the A and B sublattices, while the E\ 
representation is antisymmetric. 

Emergent 'pairing symmetry at low energy: At low en- 
ergy, the effective physics in graphene can be described 
by a relativistic dispersion near the wave vectors K± = 
(0, ±j^=)( hereafter subscript '±' always denotes the val- 
ley index). In the superconducting state of graphene, we 
also have to consider the superconducting orders near 
these vectors at low energy, which leads to the effective 
superconducting orders. In particular, when the pairing 
is between two sublattices, the effective superconducting 
orders can have new pairing symmetry around the Dirac 
cones. To see this, consider a translationally invariant 
superconducting order defined on the links of the near- 
est neighbor sites between the A and B sublattices. In 
real space, this pairing order is described by three in- 
dependent values (Ay , Ay 2 , Ay^ ) as shown in Fig[TJ In 
momentum space, the superconducting order is given by 

a 

At low energy, near the Dirac cones, the effective super- 
conducting order is given by A±(q) = A(K± +q). Given 
a small q, we have 

A ± (q) = A(K ± ) + iq- (]T 4A 4 e ±iJ? ±-H (2) 

a 

Let us consider two specific cases. The first case is ex- 
tended s-wave pairing. In this case, Ay = A, a = 1, 2, 3. 
The first term in the right side of Eq. ^ vanishes and it 
is easy to show that A±(q) = — |A(±g a + iq x ), which 
becomes a p-wave like pairing order. Therefore, the ex- 
tended s-wave pairing order in graphene at low energy is 
described by two p + ip pairing orders that are connected 
with each other by time reversal symmetry. This case has 
been studied in Ref. [H. The second case is d x 2_ y 2 + id' xy 
wave pairing on which this paper is focused. In this case, 
Ay = Ae 2w7r/3 ,a = 1,2,3. The effective superconduct- 
ing orders for small q in Eq. © become 

A d + (q) = 3Ae 1 ^ 
A d _(q) = ^A(iq x + q y ) (3) 



The first equation, Ai_(q), corresponds to s-wave pairing, 
and the second, Al (q), to p+ip-w&ve pairing. Therefore, 
at low energy, the d x i_ y i + id' xy wave pairing state in 
graphene is a superconducting state with mixed s and p + 
ip pairing orders. 

Lattice model and the quasi- particle spectrum in mean 
field: The graphene system is composed of two sublat- 
tices which are labelled as A and B, as shown in FigQ] 
If the superconducting pairing is betwen two sublattices, 
the pairing Hamiltonian can be written at the mean field 
level as follows, 

H = "* zZ^ A l B i + ^ + H.C] 

i,a,(T 

+ zZ^sSA^l^-AlBl^ + H.C] 

i,a 

- /.D4^+^W' (4) 

where the index i sums over sites on the A sublattice. 
A\ a and Bj a are creation operators for two sublattices 
and a =t)| are spin indices. The first term describes 
free band where t ~ 2.8eV is the nearest-neighbor hop- 
ping constant. In the pairing term, Ay is the spin singlet 
bond pairing order parameter which has d+id' symmetry 
under the point group Dq, i.e., Ay = Ae 2m7r / 3 , where A 
is the pairing strength. The phase of the order parameter 
winds by 4-7T around each hexagonal plaquette(shown in 
FigHJ • This ansatz preserves the rotational and transla- 
tional symmetry of the original lattice but breaks time 
reversal symmetry(TRS) manifestly. The chemical po- 
tential, /x, can be tuned by the gate voltage. 

In momentum space, we can rewrite the Hamiltonian 
in the form: H = $1 + const, where we defined 

the Nambu spinor = (A^, B^,A' ~ , B^ - ) and the 

4x4 matrix is: 

/ —fi f(k) A(fc) \ 

- = f{k)* - M A(-fc) 

k A(-fc)* fi -/(-*)* 1 ' 

\A(fc)* -f(-k) [x ) 

where the function f(k) is defined by f(k) = —t ^2 a e lk ' Sa 
and A(k) is defined by Eq.{T]). 

At low energy, we can linearize the mean field Hamil- 
tonian Eq.([ni) near the two inequivalent BZ corners K±. 
Near K±, f(K± + k) can be expanded as: 

f±(k) = f{K± + k)= v{ik x ± k v ) (6) 

where we introduced a valley dependent function f± . The 
velocity of the Dirac particles is v — y. Note hereafter 
k x , k y always refer to the relative vectors measured from 
K±. By substituting Eq.© and Eq.© into Eq.©, we 
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obtain 

H±(k) = ( v ( ±k y a * ~ k * a y) ~ V A±(fc) 

V A* _ v(±k y cr x - k x (T v ) 

where a XjV refer to the Pauli matrices. The linearized 
pairing matrices A±(fc) for the '±' valleys take the form, 

and A_(«) = A + (— /c) T . It can be easily checked that 
H+(k) and H^.{—k) transform into each other under spa- 
cial inversion, ^ Here we note that the d + id' 

pairing ansatz break TRS but preserves inversion sym- 
metry so that the valley degeneracy is unbroken. 
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FIG. 2: The energy gap -E gap as a function of chemical poten- 
tial for pairing strength A=0.001t and A=0.005t, i.e., ~3mev 
and ~15mev separately. Notably, E ga _ p is linear at low doping 
region and saturates to a constant J5|ap when fi ^> A. 

The elementary excitation spectrum can be obtained 
through Bogoliubov diagonalization. Furthermore, we 
can find the energy gap Eg ap corresponding to the mini- 
mum excitation energy. It can be shown rigorously that 
for "C A, E gap — \i. In Fig.[2]we plot the gap as a func- 
tion of chemical potential. E gap is linear in the fi <C A 
(low doping) region and saturates to a constant E g ® p for 
fx ^> A. This unique dependence of the energy gap on the 
chemical potential in the d + id' superconducting state 
stems from the mixture of the s-wave and the p + ip wave 
components. Since the p + ip wave component dominates 
at low doping, the gap depends linearly on the chemical 
potential. Since the s-wave component dominates in the 
high doping region, the gap saturates above the s-wave 
superconducting order parameter strength. 

Andreev conductance through S/N junction: In the fol- 
lowing, we show that the mixture of the s-wave and p+ip- 
wave in the d + id' superconducting state of graphene 
results in a distinctive signature in the Andreev conduc- 
tance. Consider a S/N graphene junction with the x<0 
region being the graphene d + id' superconductor and 
the x>0 region being the normal state of graphene. We 
assume that the electrostatic potential on the S side is 



lower than that on the N side by a value Uq > 0, which 
can be fixed through the gate voltage or by doping. A 
large Uq implies a heavily doped superconductor. Due to 
the spin and valley degeneracy, we can restrict the inci- 
dent state from N side to be spin up and from valley '+' 
and multiply the conductance by 4 at the end. 

Under certain voltage bias V, we expect the incidence 
of a particle excitation with energy e = eV from the x>0 
side into the junction at x=0. The general form of the 
incident wave function is ^| = <I> e (— k x , k y )e l ^~ kxX+k ^ 
where k x (k y ) is the longitudinal (transverse) component 
of the wave vector. In the scattering process, we assume 
energy and the transverse component of the wave vec- 
tor is conserved. The reflected states can be either an 
electron state ^f.=^ e (k x , k y )e l< - k:cX+k « v " > or a hole state 
i&h=& h (k' x , fc v )e^ fe « x+fc "^ where k' x is determined by 
vk'=\e — Ep\, it is negative for e < Ep (retro-reflection) 
and positive for e > Ep (specular reflection) Q • It can 
also be imaginary if vk y > |e — Ep\ and the correspond- 
ing hole state is an evanescent state near the boundary . 
<I> e and & h are 4-component spinor eigenstates of Eq.([7]) 
on N side(for which A=0) corresponding to electron and 
hole excitations, respectively. 

On the S side, we diagonalize Eq.0 and obtain the 
Bogoliubov quasi-particle states. The general form of 
the quasi-particle states on the S side is denoted by 
& s (k x , ky)e l ( k * :X+kvy " > where k x is the longitudinal compo- 
nent of the wave vector on the S side. The 4-component 
spinor Q 3 (k x ,k y ) is called electron-like (hole-like) if the 
summation of the square of absolute values of the first 
two components is larger (lesser) than that of the last 
two components. For each e and k y , we can obtain four 
quasi-particle states. Two quasi-particle states are picked 
out among four. The chosen states satisfy one of the 
following three conditions: (1) k x is real and positive 
and $ s (k x ,k y ) is hole-like; (2) k x is real and negative 
and & s (k x ,k y ) is electron-like; (3) k s x is complex and the 
imaginary part is negative. The last case corresponds to 
evanescent states near the interface. By matching the 
wave functions of both sides at the interface x=0, we can 
obtain the reflection coefficients r and ta for states ^ 
and ty^, respectively. The quantum conductance through 
the S/N junction can be calculated using the Blonder- 
Tinkham-Klapwijk formula [l7j. 

,tt/2 

G = G / (1 - \r(eV, a)\ 2 +n h \r A (eV, a)\ 2 )cosada 
Jo 

(9) 

where a = tg~ 1 (-jf L ) is the incident angle and Go = 

^-N(eV) is the ballistic conductance of the graphene 

sheet with density of states N(eV) = (£F +^ y)W (W is 
the width of the graphene sheet), nh equals 1 if the hole 
state on the N side is propagating, and it is if the state 
is evanescent. 

For ease of comparing our results with the s-wave re- 
sults in 0], we depict the normalized quantum conduc- 
tance G/Gq (as a function of bias voltage) of the S/N 
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FIG. 3: The normalized quantum conductance of a S/N 
graphene junction is shown in (a) with Ef > -E|a P an d 
(b) with Ef < Eg2p for heavily doped superconduc- 
tor(U=0.1t,i.e., 300 meV). Also shown in (c) the normal- 
ized conductance for zero bias voltage for three kinds of pair- 
ing order parameters, i.e., conventional s-wave, extended s- 
wave(bond pairing) and d + id' wave(bond pairing) 

junction with the S side being heavily doped supercon- 
ducting graphene for two cases, i.e., for Ep > E^ and 
E F < E™* in Fig. [3£a) and Fig. Mjo) respectively, where 
is the saturated gap for fi » A shown in FigO 

For Ep > Eg 1 ^, G monotonically decreases in the re- 
gion eV < A and saturates to a constant value quickly 
as eV > A. The saturation value slightly decreases with 
Ef. For Ef < -E|api the line shape is similar to the s-wave 
case, except that the unbiased conductance is nearly 2 in- 
stead of |. It is noteworthy that G is always zero at the 
point Ef — eV in Figj3j3, since there is no Andreev hole 
reflected back at this point for any angle of incidence. 

The most remarkable difference between the G/Go — 
eV curves for the conventional s-wave caseQ and the 
d + id' wave case in this paper is the value of the un- 
biased conductance, i.e., | for s-wave and nearly 2 for 
our case. In refQ, the lines are calculated in the large 
Uq limit. To make things more clear, we calculated the 
unbiased G/Gq as a function of Uq with several different 
choices of Ef and A values. In Fig.[3jc), we plot a typical 
comparison for three kinds of pairing order parameters. 
The results for the conventional s-wave(0) and extended 
s-wave ( 13]) cases show little difference. For both cases 
the unbiased G/Go quickly converges to the value of 4/3. 
But for the d + id' wave case considered here, G/Go de- 
creases slowly from 2 and converges to 4/3 after Uo > lOt, 
which is far beyond the single band edge. The most fun- 
damental difference between the d+id' pairing ansatz and 
others is that it breaks time reversal symmetry and has 



an emergent mixed s-wave and p-wave pairing at low en- 
ergy. It would be interesting to mention that in a recent 
conductance measurement on a S/N/S structure, which 
is a realization of Andreev billiards, the Andreev con- 
ductance is always peaked at zero voltage bias[16|. This 
result is consistent with our calculation and may shed 
new light on the superconducting pairing symmetry of 
the graphene. 

Realization of d + id' -wave superconducting state in 
graphene: It has been shown that the d + id'-w&ve su- 
perconducting state in graphene is a natural mean field 
solution in the presence of strong electron correlation [Toj. 
Actually, the honeycomb lattice is closely related to tri- 
angular lattice (with same lattice rotational symmetry) 
in which the superconducting state is believed to be of 
d ± irf'fllL [T3 | . Although the electron correlation in 
graphene is probably not strong enough to produce a 
d + id' superconducting by itself, it is possible to realize 
the d + id' superconducting state by including the prox- 
imity effect through a connection to another supercon- 
ductor. For instance, let us consider the geometry where 
the zigzag edge of graphene sheet is laterally connected 
to the [110] direction of a d-wave pairing superconductor 
on a square lattice (more specifically, a high Tc cuprate 
superconductor). In the coordinate system we used to 
discuss the symmetry of graphene, the order parameter 
in the dSC side is d xy . The proximity effect is therefore 
expected to induce a d xy component in the graphene side 
near the interface, which should be continually evolved 
into a d + id' pairing symmetry in the bulk due to the 
presence of electron correlation. Another possible real- 
ization is to put a dSC on top of a graphene sheet. In 
the long wave length description, the effect of lattice mis- 
match is irrelevant and d + id' -wave superconducting 
state can also be induced. More radically, it is possible 
that even a s-wave superconductor may induce d + id' or- 
der as well, as the conductance measurement [l6j that we 
mentioned earlier indicates. The immediate consequence 
of the presence of the d + id' order is the spontaneous 
supercurrent along the interface. A self-consistent study 
of such proximity effect will be presented elsewhere. 

In summary, we have shown that a combination of s- 
wave and p + ip-w&ve pairing order parameters emerges 
at low energy in the d x 2_ y 2 + id' xy spin singlet pairing 
superconducting state in graphene, which can be gener- 
ated by electronic correlation as well as induced through 
a proximity effect with a superconductor. This mixture 
of s-wave and p+ip results in distinctive superconducting 
gap functions and novel behavior of the Andreev conduc- 
tance. 
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